HOW TO KICK A SOLITON ? 

Jacek Dziarmaga* and Wojtek Zakrzewski^ 
Department of Mathematical Sciences, University of Durham, South Road, Durham, DH1 3LE, United Kingdom 

(June 25, 1997) 



A simple method how to study response of solitons in 
dissipative systems to external impulsive perturbations is de- 
veloped. Thanks to nontrivial choice of small parameter, the 
perturbative scheme captures genuine nonlinear phenomena. 
The method is developed and tested by numerical simula- 
tions for kinks in 1+1 dimensions and for skyrmions in 2+1 
dimensions. Extension to models including second order time 
derivatives is discussed. 

Let us consider a dissipative version of the <p theory 
in one spatial dimension, which is denned, in appropriate 
dimensionless units, by the nonlinear field equation 



dt<j> = dl<\> + 2[1 - 4> 2 ](f> + r)(t, x) 



(1) 



In the absence of noise, r](t, x) = 0, Eq.(|3) admits static 
kink solutions 



<j)(t,x) = F(x) = tanh(x) 
Antikinks are given by F(~x). 



(2) 



Spectrum of kink excitations. 

Small perturbations around the kink take the form 



-it 



g(x). Linearization of Eq.(|l]) with respect to g(x), 



for r](t, x) = 0, gives 

r s d2 i s ri 6 

ig{x) = --r^g{x) + 4 -j— 

oar cosh (x) 



(3) 



which looks like a stationary Schrodinger equation for a 
potential well problem. The eigenvalues and eigenstates 
can be tabulated as |Q 

0, F'(x) = ' 



3, B{x) 



cosh 2 (x) 
sinh(x) 
cosh 2 (a;) 



,2 / \ ifcrr, 3i/ctanh(x) - Stanh 2 ^^), 
A + k 2 ,g k {x) = e lkx [\ + T^tf — 1 ' ( 4 ) 

where k is a real momentum. The zero mode (7 = 0) is 
separated by a gap from the first excited state (7 = 3). 
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Naive perturbation theory (NPT). 

Let us consider the response of the kink to a soli- 
tary perturbation i](t,x) — 5(t)ip(x). Let the unper- 
turbed state be <f)(t < 0, x) = F(x)\ integration of Eq.(Q) 
over time in an infinitezimal neighborhood of t — re- 
sults in the initial configuration (f>(0 + , x) = F(x) + ip(x). 
The function tjj(x) can be expanded in the basis (^) as 
= $,F'(x) + 4> B B(x) + dk ${k)g k (x). If yj(x) 
is sufficiently small, such that Eq.(0) can be linearized 
around <fi(t, x) — F(x) for t > 0, then the perturbed field 
is well approximated by 

4>(t > 0, x) w F(x) + i> z F'(x) + e- 3t 4> B B(x) + 

dk e- (4+k2 ^yl(k)g k (x) . (5) 



After a sufficiently long time, as compared to the relax- 
ation time of the first excited state equal to 1/3, the field 
relaxes to 



oo, a;) w F(x) +ip z F'(x) 
_ Jdx F'(x)ip{x) 
■ ~ Jdy F'(y)F'(y) ' 



F[x + vl z 



(6) 



The kink is shifted to a new position, the shift is given 
by a projection of the initial perturbation tp{x) onto the 
kink's translational zero mode. However plausible it may 
sound, this result is not very accurate even for small per- 
turbations. The reason is that the spectrum (^) is rel- 
evant for perturbations around the initial kink at the 
origin but it becomes less relevant as the kink drifts to 
its final location at — ip z . The excited states in (|J) are 
not orthogonal to the zero mode F'\x + tp z ] of the final 
kink. 



Improved perturbation theory (IPT). 

An appropriate choice of perturbation theory may im- 
prove the accuracy. Let us consider again the single kick 
rj(t, x) = 5(t)tp(x) and let us rewrite the initial perturbed 
field in the following form 

<£(0 + , z) =F(x)+i/j(x) =F{x-£+)+ti{0 + ,x) . (7) 

The above equation is just a tautology, unless we impose 
a constraint on the function /j,(0 + ,x). We require the 
function [i(0 + ,x) to be orthogonal to the zero mode of 
the kink at £ + , 



1 



= 



dx F'(x-£ + )n(0 + ,x) 



(8) 



This orthogonality implies that /i(0 + , x) can be expressed 
as a combination of the excited states in the complete set 
(@) shifted to £+, 

/+oo 
dfc/i(fc). 9fc (x-e + ) • 
-oo 

(9) 

So far we have made no approximation. The perturbative 
approximation is defined by the assumption that /i(0 + , x) 
is sufficiently small, so that Eq.(Q) can be linearized in 
/i(t, x) around F(x — The linearized solution is given 
by 

H(t >0,x)= fi B e- 3t B(x - £+) + 

dkfi(k)e- {i+k2)t 9k(x-e) ■ (10) 



After a sufficiently long time fJ,(oo, x) = and the initial 
condition on the R.H.S. of Eq.(^) evolves into F(x — 
Thus £ + is the final position of the kink after perturba- 
tion. 

To find £+, let us project Eq.(@) on the zero mode 
- £+). Using Eq.(|) and the fact that F(x) is or- 
thogonal to F'(x), we obtain 



-(-oo 



dxF'(x-^ + )[F(a;)-|-V'(a;)] =0 , (11) 



which can be further rewritten as 

r° dI (12) 

J-oo cosh (xj 

where U is the smooth, antisymmetric and monotonic 
function 



u(e) = 2[ 



i 



tanh(f) sinh 2 (£) 



(13) 



Already from the general form of (|T^) it is clear, that 
£ + is not linearly dependent on the magnitude of the 
perturbation ^>{x). This response is nonlinear although 
we use a sort of perturbation theory. 



Numerical simulation. 

Before we proceed with general discussion, let us con- 
sider some quantitative examples and compare IPT pre- 
dictions with direct numerical simulations of Eq. (jl|) . We 
take the perturbation in the form 



if) (a 



cosh [x — y) 



(14) 



where A is the amplitude and y the localization of the 
perturbation. Eq.(|12|) takes the form 



A = 



sinh 3 (£+-y)[ 



tanh(f"t 



2[{t+-y) cosh(r 



V) 



sinh(£+ - y)] 



(15) 



Fig.l illustrates the dependence of the response £ + on 
the amplitude A for a central perturbation with y = 0. 
IPT works well up to A 1 , or up to soliton shifts compa- 
rable with soliton size. Fig. 2 shows the data for A = 0.2 
and y varying in the range —2 < y < 2. If the perturba- 
tion hits to the right of the kink, y > 0, and is such that 
the kink is driven to the right, A > 0, the response is 
amplified because the perturbed kink is flowing towards 
the source of the perturbation. This positive feedback, 
which is missing in NPT, explains the asymetry in Fig. 2. 

For some y Eq.(|l^) does not have a unique solution 
£ + = £ + (A) in the whole range of A. This usually hap- 
pens for \A\ > 1, when the perturbation theory is not 
reliable a priori. As an example, Fig. 3 illustrates the re- 
sponse for y = 2 and A > 0. The theoretical function 
is not unique around A ps 1.3. Numerical results coarse 
grain over this singularity. The perturbation (|l4j) with 
y = 2 creates a virtual antikink-kink pair just to the right 
of the original kink (^|). The result is a kink- antikink-kink 
configuration. For a weak perturbation (A < 1.3) the an- 
tikink annihilates with the right kink, leaving the original 
left kink shifted a bit to the right. For A > 1.3 the an- 
tikink annihilates the original left kink, the net result is 
a huge shift to the right (£ + > A) , much larger than that 
expected from NPT. This is again the positive feedback 
of Fig. 2 but in a caricaturely enhanced form. This is a 
genuine nonlinear response and it is amusing to find that 
IPT does capture its essential features. 



Response to subsequent perturbations. 

£ + is the position the kink relaxes to after a sufficiently 
long time. Before the final equilibrium is reached, the 
function /i(t,x), see Eq. (10), carries information about 
the initial perturbation. This information is lost after 
roughly 1/3 units of time, which is the relaxation time 
of the breather B(x) (||). As an example, let us consider 
the kink perturbed by the noise 



r)(t, x) = S(t)ip(x) + S(t - T)^{x) 



(16) 



where *p(x) is given by ( |l4| ) with, say, y — and A > 0. 
At the time t = T the field is hit for the second time, 



4>(T+, x) = F(x - £+) + M (T, x) + t(x) 
F(x-e + )+ft(T + ,x) , 



(17) 



while it still remembers the first perturbation /it(T, x) 
Repeating steps leading to (|l2|), one obtains 



2 



cosh (x) 



(18) 



There are two limiting cases. If T = the two kicks 
add up to a solitary perturbation 2?/>(x). If the function 
£, + (A 7 y) is a solution of Eq. (|l5|) , then the final shift of 
the kink can be expressed as £ ++ (T = 0) = ^+(2A, 0). 
In the limit T = oo, the kink has enough time to set- 
tle down at £ + (A, 0) after the first perturbation before 
it is shifted by the second kick to £ ++ (T = oo) = 
£+(A,0) + f+[A, -£+(A,0)]. The two shifts are differ- 
ent, £ ++ (T = oo) > £(T = 0). Fig. 4 gives numerical 
results for various values of T. The limit £ ++ (T = oo) 
is essentially achieved for T greater than the relaxation 
time of the breather mode ([|), T > 1/3. For shorter in- 
tervals T the memory function /i(T, x) in Eq. (|l8|) cannot 
be neglected. 

The second kick probes the position of the kink at the 
time T after the first perturbation. Roughly speaking, 
the kink remains unshifted for T < 1/3, it jumps to 
£+(A,0) for T > 1/3. The <5(t)-like time shape of our 
impulse may be good approximation for impulses lasting 
less than 1/3. 



Higher order time derivatives. 

IPT can be generalized to models with higher order 
time derivatives. For example, a modification 



dtfi — > d t cj> + adfc/) 



(19) 



on the LHS of Eq. (|T|) changes the initial conditions af- 
ter the perturbation r](t,x) = 6(t)i/}(x) to: </>(0 + ,x) = 
F(x) , <9 t 0(O + ,x) = tp(x). The perturbed field can be 
written as <p(t > 0,x) = F(x~£,+) + (A 1 +A 2 e- t / a )F l {x- 
£ + ) + R(t,x). R(t,x) can be expanded in the excited 
states which decay with time. For £ + to be the final kink 
position, <f>(po,x) — F(x — £ ), we must demand A\ = 0. 
This constraint plus the initial conditions result in the 
following equation for £ + 



'dx^ + r) =0 
cosh (x) 



(20) 



which is very similar to Eq.( [T2|). It is important to realize 
that the a — > limit of Eq.(]20|) docs not give Eq.(12), 
which holds for a = 0. The more general model is dis- 
continuous at a = 0. 

We have performed several numerical simulations to 
test Eq. (|20|) . The accuracy is the same as for a = 0. 



IPT and skyrmions in 2+1 dimensions. 

Once IPT has passed the tests for kinks, it would be 
interesting to check how it works in higher dimensional 



models. Let us consider the diffusion equation in the 
x = (xi, X2) plane 



5W 
SM 



d t M = P £l [-^-+f j (t,x)] 



(21) 



where M(t, x) is a unit (MM = 1) magnetization vector, 
P A j is a projector on the subspace orthogonal to M, W — 
J d 2 x w is the static energy functional with the energy 
density 



w = ^d k MO k M + ^ ai (Ml + Mi) + 



1 



a 2 [(d k Md k M) 2 - [d k MdiMf 



(22) 



being a sum of the exchange, easy axis anisotropy and 
Skyrme energy densities respectively. The energy is min- 
imized by the skyrmion configuration [p[ 

F(x\,X2) — {sm[q(r)} cos#, sin[g(r)] sin0, cos[q(r)]} (23) 

where (r, 9) are polar coordinates in the (xi,X2) plane 
and q(r) is a profile such that q(Q) — ir and q(oo) = 0. 
q(r) has been determined by a shooting method for a\ 
I and a 2 = 0.1, the result is shown in Fig. 5. 

The initial field after the perturbation fj(t,x) = 
S(t)ip(x) is a bit complicated because of the projector 
on the RHS of Eq. (pl|) . At any point x the integration 
of Eq.(|2l]) around t = gives the initial perturbed mag- 
netization 




(24) 



Similarly as for kinks, compare Eqs.(|7|,|ll|), the final po- 
sition £ + of the skyrmion can be found from the condi- 
tion that the initial perturbed field is orthogonal to the 
k = 1, 2 transnational zero modes of the final skyrmion 



= / d z x M(0^ 



(25) 



We performed numerical simulations of Eq. (pl|) for 
the kick V>(x) = "O-Sgf^-^Xx* — y) with y — (y%, 0) re- 
stricted to the xi-axis. As a result of the perturbation the 
skyrmion was shifted to a new position at £ + = (£i~,0). 
The function g£ (yi) is presented in Fig. 6. The maximum 
of the function is shifted to the right and the minima are 
shifted to the left as a result of similar positive feedback 
as in Fig. 2. 

IPT works well for skyrmions in 2 + 1 dimensions. The 
model @ can be generalized to, say, Landau-Lifshitz 
equation with dissipation. Modification of the LHS of 
Eq.©, d t M d t M + AM x d t M, changes M(0+,x) 
but the condition (pa) remains the same. 



3 



Remarks. 



The theory outlined above can be upgraded to a fully 
fledged theory of soliton diffusion driven by a random 
noise. Another possible line of development is a general- 
ization to the case of dissipative multisolitons H . 
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FIGURE 1. for the kink as a function of A for 
y = 0. Dashed line - NPT, solid line - IPT and dots - 
numerical simulation. 
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FIGURE 3. £ + for the kink as a function of A for FIGURE 4. Numerical simulation results for £++ af- 
= 2. ter two subsequent perturbations separated by T as a 

function of T. Amplitude A=0.3. 
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FIGURE 5. The magnetization M 3 = cos[g(r)] as a 
function of r. 
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FIGURE 6. £i~ for the skyrmion as a function of y\ for 
A = 0.3. Dashed line - NPT, solid line - IPT and dots - 
numerical simulation. 
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